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ABSTRACT 

Context. Establishing the multiplicity of O-type stars is the first step towards accurately determining their stellar parameters. 
Moreover, the distribution of the orbital parameters provides observational clues to the way that O-type stars form and to the in- 
teractions during their evolution. 

Aims. Our objective is to constrain the multiplicity of a sample of O-type stars belonging to poorly investigated OB associations in the 
Cygnus complex and for the first time to provide orbital parameters for binaries identified in our sample. Such information is relevant 
to addressing the issue of the binarity in the context of O-type star formation scenarios. 

Methods. We performed a long-term spectroscopic survey of nineteen O-type stars. We searched for radial velocity variations to 
unveil binaries on timescales from a few days up to a few years, on the basis of a large set of optical spectra. 

Results. We confirm the binarity for four objects: HD 193443, HD 228989, HD 229234 and HD 194649. We derive for the first time 
the orbital solutions of three systems, and we confirm the values of the fourth, showing that these four systems all have orbital periods 
shorter than 10 days. Besides these results, we also detect several objects that show non-periodic line profile variations in some of 
their spectral lines. These variations mainly occur in the spectral lines, that are generally affected by the stellar wind and are not likely 
to be related to binarity. 

Conclusions. The minimal binaiy fraction in our sample is estimated to be 21%, but it varies from one OB association to the next. 
Indeed, 3 O stars of our sample out of 9 (33%) belonging to Cyg OB 1 are binary systems, 0% (0 out of 4) in Cyg 0B3, 0% (0 out of 
3) in CygOB8, and 33% (1 out of 3) in CygOB9. Our spectroscopic investigation also stresses the absence of long-period systems 
among the stars in our sample. This result contrasts with the case of the O-type stellar population in NGC 2244 among which no 
object showed radial velocity variations on short timescales. However, we show that it is probably an effect of the sample and that this 
difference does not a priori suggest a somewhat different star forming process in these two environments. 

Key words. Stars: early-type - Stars: binaries: spectroscopic - Open clusters and associations: individual: Cygnus OB 1 - Open clusters 
and associations: individual: Cygnus 0B3 - Open clusters and associations: individual: Cygnus OB8 - Open clusters and associations: 
individual: Cygnus OB9 
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CO 1. Introduction renders our knowledge of massive stars still fragmentary. In 

^, , . ,. . ^ . . ^ , addition, the distribution of orbital parameters allows us to 

- The multiphcuy of massive stars constitutes one of the most es- -^^- ^^ -^^^ ^^^^^ aveniom such as whether the binaiy 

:^ ' sential ingredients for understanding these objects Establishing fac tion of m assive st ars can be related to the stellar density 

; itis the first step m obtaimng valuable physical parameters, (Pen l^7ir^ ri993: Gar cia & Mermilli"^ ^00 11. The orbital pa- 

. such as their masses, radii and luminosities. The binanty also ^^^^^^^^ ^^^^- ^^ ^ ^^^^^ ^^^^ ^^^-^^ important clues 

_ C3_ affects the evolutionai-y paths of the components in such sys- conditions prevaiUng in the birthplace during their for- 

tems in comparison to single stars through phenoin ena, such a s ^^^.^^ dynamical interations that occur during the 

tidal interations or Roche lobe overflows (Sana & Evans|M). eariiest stages of their evolution. Several studies have already 

Knowing this multiplicity also allows us to be more accurate on ^^^^ performed in spectroscopy to determine the binai-y frac- 

the wind parameters, especia ly the mass-loss rates, or on the in- tion of O-type stars in nearby young open clusters (see e.g., 

trinsic X-ray luminosities. Moreover, the distribution of the or- |Sana &Evansl l20Tll and references therein). All these anal- 

bital parameters for the massive systems provides us with ad- y^es revealed a n averaged binary fraction of about 44 ± 5% 

ditional information on the way that massive close binaries ai-e jsana & EvansllIoTlh . However, to investigate such a large pa- 

■ . r , r ■ r , , Tameter space. It is necessai^ to comblne Other obseFvatlon tech- 

Many different scenarios for the formation of these ob- ^. ^^^^ interf erometry ([N elan et al. 20 Q4|), speckle in- 



jects have been proposed (see Zinnecker & Yorke 2007, for a terf^erom etrv (see e.g.. ^ason et al 1998: Maiz-Apellaniz et al 

complete review), but none of them totally explains what hap- 2004; M ason et alJ 120091) . or adaptative optics dTurneret" 



pens during their formati on or during its eariiest stages. This 20O8). Indeed, all these large surveys give the advantages of bet- 
Send ojfprint requests to: L. Mahy ^r statistics whilst the surveys devoted to young open clusters 
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or OB associations allow mor e homogeneous stellar populati ons 
to be investigated as noted bv lKkninki & KobulnickM j2012l) . 

The Cygnus area is an active star-formin g region in the 
Milk y Way. Its relative proximity (0.7 - 2.5 kpc. lUvaniker et aU 
I2OOI ) makes this region suitable for studying stellar populations 
and star formation processes (see Reipurth & Schneider 2008, 
for a review). The Cygnus constellation harbours nine OB as- 
sociations and at least a dozen of young open clusters. To bet- 
ter constrain the environment of the Cygnus region. Fig JT] rep- 
resents a schematic view with the locations, in galactic coor- 
dinates, of the different OB associations and of several young 
open clusters. The central association Cyg OB2 is ceitainly the 
most famous and one of the youngest of the Cygnus region 
with CygOBS. To be more accurate, CygOB2 has in fact been 
shown to possibly have two populations with ages of 2 - 3.5 and 
5 Myrs (Wright et al. 2010), whilst the O stars seem to belong 
to a y ounger population, possibly aged about 2 Myrs (Hanson 
l2003h . In its suiToundings, the associations Cyg OBI, CygOB3, 
CygOB7 and CygOB9 are old er with ages of about 7.5, 8.3, 
13.0 and 8.0 Myrs, respectively (lUvaniker et al.ll200 r). Though 
the Cygnus region is rich in O-type stars, its heavy absorption 
prevents one from obtaining a better estimate of the stellar popu- 
lation even though about 100 hot stars are expected i n Cyg OB 1 , 
CygOB8 and CygOB9 together (iMel'N ik & Efre movl Il995h 
and about 10 O- stars are expec ted only i n the CygOB2 asso- 
ciation (Kno dlseder et a l. 2002; Comeron e t al. 2002). However, 
the membership of certain massive stars to a given cluster or as- 
sociation is sometimes not clearly defined due to the complexity 
of this area. Recent investigation of the multiplicity of massive 
stars in CygOB2 association peiformed by Kiminki et al. (2012, 
and subsequent papers) has shown that the hard minimum binary 
fraction for the massive stars in this association is estimated at 
21%. Such an investigation for the massive star population in 
the other OB associations in the Cygnus complex is currently 
lacking. 
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Fig. 1. Schematic view of the different OB associations and 
young open clusters in the Cygnus region. 



We have undertaken a long-term spectroscopic survey of 
nineteen O-type stars sampled over four different OB associa- 
tions (Cyg OBI, CygOB3, CygOB8 and CygOB9) as well as 
in several young open clusters belonging to these associations 
(e.g., Berkeley 86, NGC6871, NGC6913,...). Their membership 
to a given as sociation or c luster is established on the basis of the 
catalogue of lHumphrevsl (l978). However, many other O-type 
stars belong to these associations but we only focus on these 



targets because they are relatively bright and their binary status 
was not yet subject to a detailed spectroscopic study, except for 
HD 229234. We will also include other stars already known as 
binaries in the discussions on the binary fraction in the different 
Cygnus OB associations. Moreover, the majority of the objects 
investigated in the present paper were poorly observed in the past 
mainly due to the high extinction towards the Cygnus region. In 
this context, this paper (the first of a series of two) aims at es- 
tablishing the multiplicity as well as the spectral classification of 
these stars. This analysis focuses on other region s in the Cygnus 
complex than the large survey of iKimrnki et al.l (1201 2l and the 
subsequent papers) which was devoted to the Cyg OB2 associa- 
tion. Therefore, this analysis can thus be seen as a complement 
to the study of the binary fraction in Cygnus region. 

The present paper is organized as follows. The observing 
campaign is described in Sect.|2] The binary status of the stars in 
our sample is discussed in the next four sections according to the 
different OB associations. Sectionals devoted to O-type stars in 
the Cyg OBI association. Sect. |4] to these objects in CygOB3, 
Sect.|5]to CygOB8, and Sect.|6]to CygOB9. Finally, Sect.Qdis- 
cusses the observational biases and the binary fraction that we 
determine whilst Sect. [8] summarizes our results. 

2. Observations and binary criteria 

Deriving the binary fraction of a given population requires an 
intense survey of the stars but also a good time sampling to con- 
strain the short- as well as the long-period systems. This moni- 
toring has to be as homogeneous as possible because combining 
data from the literature may provide eiToneous results. Indeed, 
the literature (especially the older papers) often does not specify 
on which spectral line the radial velocities (RVs) were measured. 
Therefore, large differences can be observed between RVs mea- 
sured on two different lines of a presumably single star. This 
can lead to some overdeterminatio ns of the spectroscopic binary 
fraction as notably observed in iGar cfa & Mermilli odI (j20(31 1 
for the young open cl t ister IC 1805 (Ulauw & De Becker 2004; 
iDe Becker et al.ll2006t iHillwig et aD l2006). Moreover, it is im- 
portant to keep in mind the possible observational biases that 
can prevent the detection of a binary system. Therefore, stars for 
which no evidence of binarity is found can never be definitely 
considered as single. Even if no RV shift is detected, the system 
could be seen under a paiticular orientation, have a very long 
period or perhaps a high eccentiicity, thereby making the RV 
variations not significant over a long timescale. 

We performed a spectroscopic survey spread over three years 
(from 2008 to 2011) of all the nineteen stars of our sample, 
except for HD 194280 which was observed between 2004 and 
2007 (see Table[TJ. This selection of stars however constitutes 
another observational bias because we limited our observations 
to stars bright enough to yield good quality data with exposure 
times not longer than one hour with 1.5 or 2- meter class tele- 
scopes. Considering the high - and heterogeneous - extinction 
towards the Cygnus region, the present monitoiing of O-type 
stars (among which some could also be member of multiple sys- 
tems) is far from being complete. Neveitheless, we collected a 
set of 274 spectra for these nineteen stars by using three dif- 
ferent instruments. Table [1] lists the stars and the parameters of 
the observing runs. Table |2] available electronically, gives the 
Heliocentric Julian Dates (HJD) and the RVs measured on dif- 
ferent spectral lines for all the observations. For the SB2 sys- 
tems, the RVs refined by cross-correlation are also provided. In 
the latter table as in the following of the present paper, HJD is 
expressed in format HJD - 2 450 000 for readability. Our dataset 
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thus allowed us to study the RV and line profile vaiiations in the 
observed spectra over short- (a few days) and long-term (a few 
years) timescales. 

A first part of our dataset was obtained with the Aurelie 
spectrograph at the 1.52m telescope at Observatoire de Haute- 
Provence (OHP, France). The spectra are obtained with a resolv- 
ing power of about R - 9 000, they are centred on 4650 A and 
cover wavelengths between 4450 A and 4900 A. However, for 
some spectra of HD 193443, another configuration of the instru- 
ment was chosen. In this case, the resolving power was equal 
to R = 12000, these spectra are centred on 4550 A and cover 
4450 to 4650 A. Typical exposure times between 25 and 60 min 
yielded spectra with a signal-to-noise ratio larger than 150 as 
measured in the continuum close to 4800 A. The data reduction 
procedure, performed with the MIDAS software, is described in 
iRauw &De Beckeil (12004 . 

Another part of the data was collected with the ESPRESSO 
spectrograph mounted on the 2.12m telescope at Observatorio 
Astronomico Nacional of San Pedro Martir (SPM) in Mexico. 
This echelle spectrograph provides spectra spread over 27 or- 
ders, covering a wavelength domain between 3780 A and 6950 A 
with a spectral resolving power of R = 18 000. Typical exposure 
times ranged from 15 to 30 min. We combined the consecutive 
exposures of an object in a given night to obtain signal-to-noise 
ratios close to 150. As a result, our time series did not allow us 
to investigate RV variations with timescales shorter than 0.8 day, 
which is not a problem for multiplicity investigations as typically 
the shortest periods for massive stars are longer than one day 
(as Kepler's third law would imply orbital separations smaller 
than the radii of the components for such short periods). The 
data were reduced using the echeUe package available within the 
MIDAS software. 

We also retrieved eight spectra taken before 2008 from the 
Elodie archives to complete our dataset. Elodie was an echelle 
spectrograph, mounted at the 1.93m telescope at OHP, with 67 
orders in the [3850 - 6850] A spectral domain and a resolving 
power of 7? = 42 000. 

We fit one or two Gaussian profiles on several strong spec- 
tral lines to estimate the RVs of the stars when the line profiles 
were not too much broadened by rotation. For the rapid rotators 
(v sin i > 300 km s"'), we used a synthetic line profile broad- 
ened with the projected rotational velocity of the star to measure 
the Doppler shifts by least square. This RV analysis is performed 
by using the rest wavelengths quoted by Conti et al. (1977) for 
the [4000 - 5000] A wavelength domain and UnderhilU 1994.) 
for wavelengths above 5000 A. The common wavelength region 
between the majority of the data is [4450-4900] A. In this spec- 
tral band, some interstellar lines exist but they are not sufficiently 
strong and narrow to use them as indicators of the RV unceitain- 
ties. Our previous experience with the same spectrographs (see 
iMahv et al. 2009) showed that a standard deviation on the RVs 
of about 7-8 km s"' is commonly achieved with these spec- 
tra. We therefore consider the RV changes as significant when 
the standard deviations computed from the RV measurements 
reach at least this value. For rapid rotators, we consider them 
as variable if standard deviations of at least 15 km s"' are de- 
termined. We stress that the thre shold value for the rapid rota- 
tors has been increased relative to lMahv et al.l(l2009l) because the 
present data are slightly noisier. As a consequence, we consider, 
in the present paper, that a star is an SBl binary system if the RV 
vaiiations are significant (i.e., the standard deviation is at least 
larger than 7-8 km s ' for the slow and mid rotators and at least 
15 kms ' for rapid rotators) and periodic. We are dealing with an 



SB2 if the signature of a companion whose spectral lines move 
in anti-phase with those of the main star is observed. Moreover, 
if in the cases where significant variations of the RVs are mea- 
sured on all the spectral lines without any periodic motion, the 
object will be labelled as a "binary candidate". Finally, we also 
calculate the Temporal Variance Spectrum (TVS. lFuUerton et aLl 
1996) to search for line profile variations. The TVS provides 
quantitative information of the temporal variability as a function 
of the wavelength. If the variations exceed the threshold we can 
reject the null hypothesis of non-variability of the spectral lines. 
When such variations are detected for all the line profiles, the 
TVS may then be considered as an additional indicator of multi- 
plicity and, in this case, the star is reported as binary candidate. 

The SB2 binary systems detected in the present work are 
analysed as follows. After having measured the RVs of both 
components, we use these values as an input of our disentangling 
programme. This programme is based on the lGonzalez & Levatol 
(2006.) technique and it allows us to compute the individual 
spectra of both components, and to refine the RVs by apply- 
ing a cross-correlation technique. It is thus possible to obtain 
more reliable RVs at phases where the spectral lines are heavily 
blended. The cross-correlation windows used for the RV deter- 
minations generally gather the strongest helium lines (He 1 447 1 , 
Hen 4542, Hen 4686 and Hei 4713) for both components. 
With these refined RVs, we then apply the Heck -Manfroid- 
Mersch method (h ereafter HMM, iHeck et all 1 19851 revised by 
IGosset et al.ll2001h to the time series of RVs - RVp to determine 
the orbital period of the system. The 1-cr error-bar on the pe- 
riod is estimated by assuming an uncertainty on the frequency 
of the peak associated to the orbital period equal to 10 % of the 
natural width of the peaks in the Fourier power spectrum (e.g., 
Av ~ (10 Artot)"', where Artot is the time elapsed between the 
first and the last observations of our campaign). The orbital pe- 
riod is then used to determine the orbital solution of the system 
by applying the Liege Orbital Solution Package (LOSfQ)- To fin- 
ish the analysis of the systems, we estimate, from the spectral 
classification of each component, the spectroscopic brightness 
ratio of the system. For that purpose, we compute the mean ratios 
between the observed equivalent widths (EWs) and the "canon- 
ical" EWs for stars with same spectral classifications as those 
of the compon ents of binary systems. The se " canonical" va lues 
are taken from IConti & Alschuleii (|1971|) and IContil (Il973l) . In 
addition, we associate these values to EWs measured from syn- 
thetic spectra of stars having also similar spectral classifications. 
The 1-cr error-bar given on the brightness ratio corresponds to 
the dispersions of the ratios measured on each spectral line. This 
brightness ratio is essential to correct the disentangled spectra in 
order that they are comparable to those of single stars. 

Finally, for all the O-type stars of our sample, we also de- 
ri ve their spectral type. For t h at, we use the qua ntitative c riteria 
of ' Conti & Alschuleii (fT97TI) . ICoiitil (fT97l and iMaflivsl (fT988l 
1989r These criteria rely on the EW ratio of diagnostic lines. 
These EWs are measured on the highest signal-to-noise spec- 
tra of each objects. We adopt, in the present paper, the follow- 
ing usual notations: log W - log(£'W447i) - log(£'W4542) for 
the spectral type, logW" = log(£'W4089) - log(£'W4i43) and 
log W" = log(£'W4388) + log(£'W4686) for the luminosity class. 
However, the log W" criterion is restricted to 07 - 09.7 stars 



' This program, maintained by H. Sana, is available at 
http://www.science.uva.nl/~hsana/losp.htmI and i s ba sed on the 
generalization of the SBl method of Wolfe et"an lll967h t o the SB2 
case a long the lines described in iRauw et al.l ( |200(]|) and ISana et aT] 
j2006h . 
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whilst the log W can only be used for O 8 - O 9.7 stars. To 
be applicable, this last criterion requires that the stars be single 
or that the brightness ratio between both components be known. 
The uncertainties on the spectral type determination are gener- 
ally of about one temperature and luminosity class. To complete 
this classificati on, we also use the "f tag" notation presented by 
IWalbornI (119711) and the subsequent papers. 
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Table 1. List of stars investigated in this study. Ar corresponds to the time elapsed between the first and the last observation of an observing run 
and is expressed in days, and N represents the number of spectra for each star. 



Stars 


Cyg 


Obs. run 


Instrument 


Wavelength domain 




N 


HD 193443 


CygOBl 


Aug. 2004 


OHP-Elodie/1.93m 


[3850 


- 6850]A 


- 


1 






Aug. 2005 


OHP - Elodie/1.93m 


[3850 


- 6850]A 


- 


1 






Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


7.06 


4 






Oct. 2008 


OHP- Aur61ie/1.52m 


[4450 


- 4900]A 


3.17 


3 






Dec. 2009 


OHP - Aur61ie/1.52m 


[4450 


- 4650]A 


2.98 


3 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


19.98 


13 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


5.94 


7 






Dec. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4650]A 


4.00 


5 


HD 193514 


Cyg OBI 


Aug. 2004 


OHP-Elodie/1.93m 


[3850 


- 6850]A 


- 


1 






Aug. 2005 


OHP-Elodie/1.93m 


[3850 


- 6850]A 


- 


1 






Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


7.17 


4 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.07 


5 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Jun. 2010 


OHP- Aur61ie/1.52m 


[4450 


- 4900]A 


5.83 


6 






Aug. 2010 


OHP - Aurehe/1.52m 


[4450 


- 4900]A 


4.93 


4 


HD 193595 


CygOBl 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


7.13 


4 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.14 


2 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.00 


2 


HD 193682 


CygOBl 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


7.02 


4 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.06 


2 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


- 


1 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


- 


1 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


2.98 


2 


HD 194094 


CygOBl 


Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.93 


2 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.94 


2 






Sep. 2011 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


- 


1 


HD 194280 


CygOBl 


Oct. 2004 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.05 


3 






Jun. 2005 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


8.03 


7 






Sep. 2006 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.07 


8 






Oct. 2006 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


- 


1 






Nov. 2007 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


26.95 


23 


HD 228841 


CygOBl 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


6.08 


3 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


- 


1 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


- 


1 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


- 


1 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


3.18 


2 


HD 228989 


CygOBl 


Aug. 2009 


SPM - Espresso/2.12m 


[3780 


B 

- 6950]A 


0.80 


2 






Jul. 2010 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


- 


1 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


5.13 


5 






Jun. 2011 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.94 


5 






Sep. 2011 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


1.98 


5 


HD 229234 


CygOBl 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.06 


2 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.82 


2 






Jun. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


- 


1 






Aug. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


0.96 


2 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


6.06 


4 






Jul. 2010 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


1.06 


2 






Aug. 2010 


OHP - Aurelie/1.52m 


[4450 


- 4900]A 


5.86 


6 


HD 190864 


CygOB3 


Aug. 2001 


OHP-Elodie/1.93m 


[3850 


- 6850]A 


- 


1 






Sep. 2008 


OHP - Aurelie/1.52m 


[4450 


- 4900]A 


7.05 


5 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


3.25 


3 






Jun. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


5.91 


2 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.00 


2 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.00 


4 


HD 227018 


CygOB3 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


7.05 


4 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.02 


4 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


2.02 


2 


HD 227245 


CygOB3 


Aug. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


0.86 


2 
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Table 1. continued. 



Stars 


Cyg 


Obs. run 


Instrument 


Wavelen 


gth domain 


AT 


N 






Jun. 2011 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.81 


3 


HD 227757 


CygOB3 


Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


4.98 


2 






Jul. 2010 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 




1 






Jun. 201 1 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.94 


3 


HD 191423 


Cyg OB 8 


Aug. 2004 


OHP - Elodie/1.93m 


[3850 


- 6850]A 


— 


1 






Sep. 2008 


OHP - Aur61ie/1.52m 


[4450 


- 4900]A 


7.09 


4 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


— 


1 






Jun. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


3.17 


2 


HD 191978 


Cyg OB 8 


Aug. 2001 


OHP-Elodie/1.93m 


[3850 


- 6850]A 


1.05 


2 






Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


5.94 


3 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


— 


1 






Jun. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


8.92 


2 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


4.11 


3 


HD 193117 


Cyg OB 8 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


6.82 


3 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


— 


1 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


5.05 


2 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


3.20 


2 


HD 194334 


CygOB9 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


6.75 


2 






Oct. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


— 


1 






Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Aug. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


0.87 


2 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


1.07 


2 






Jul. 2010 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


1.03 


2 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


3.87 


1 


HD 194649 


CygOB9 


Jun. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Aug. 2009 


SPM -Espresso/2. 12m 


[3780 


- 6950]A 


— 


1 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


5.95 


3 






Aug. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


2.13 


2 






Jun. 2011 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.93 


5 






Sep. 2011 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


3.98 


5 


HD 195213 


CygOB9 


Sep. 2008 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 


6.96 


3 






Oct. 2008 


OHP - Aurelie/1.52m 


[4450 


- 4900]A 




1 






Aug. 2009 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 


0.96 


2 






Jun. 2010 


OHP- Aurelie/1.52m 


[4450 


- 4900]A 




1 






Jul. 2010 


SPM - Espresso/2. 12m 


[3780 


- 6950]A 




1 






Aug. 2010 


OHP- AureUe/1.52m 


[4450 


- 4900]A 


4.91 


2 
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3. 0-type stars in the Cyg 0B1 association 

3.1. Gravitationally bound systems 
3.1.1. HD 193443 

Reported as binary bv lMuUerl (Il954l) . this star has never been the 
target of an intense spectroscopic monitoring. Therefore, no or- 
bital solution is known for the system nor are there spectral clas- 
sifications for the components. We acquired thirty-seven spectra 
of HD 193443 from Aug 2004 (HJD = 3246.4777) to Dec 2010 
(HJD = 5543.2387). These data reveal a periodic variation of the 
line profiles, clearly visible in the He 1 4471 and He 1 4713 lines. 
Even though the signature of the secondary is detected, the small 
RV separation between both components prevents us from com- 
pletely resolving the individual line profiles throughout the orbit. 
The secondary's lines are thus never totally separated from those 
of the primary even for the helium lines or the metallic hnes as 
it is shown in the upper panel of Fig. |2] 

The analysis of the SB2 system is performed as explained in 
Section |2l However, the Gaussian fit is made on the bottom of 
the absorption lines of the primary component. At the maximum 
of the RV separation, we also fit a second profile to determine 
the RVs of the secondary. The HMM technique measured on the 
RVs refined by the disentangling programme yields a period of 
7.467 + 0.003 days. From this period, we compute the orbital pa- 
rameters of HD 193443 listed in Tableland the best fit of the RV 
curves, shown in Fig.|2] is obtained for a non-zero eccentricity. 

The study of the disentangled spectra of HD 193443 shows 
that the He n 4542 line is present but is weak for both compo- 
nents relative to the He 1 447 1 line. This indicates two late O-type 
stars. We compute, from the observed spectra, log W - 0.321 
and log W - 0.483 for the primary and the secondary, respec- 
tively, thereby corresponding to respective subtypes of O 9 and 
O 9.5. However, our dataset does not allow us to constrain the lu- 
minosity classes of these stars. Indeed, the majority of our data 
does not cover the Si iv 4089 and He 1 4143 lines. Moreover, for 
those that do, the RV separation is insufficient to resolve with 
accuracy the signature of the secondary component and thus to 
correctly fit it with Gaussian profiles. To determine the bright- 
ness ratio, we measure the EWs of the He: 4471, 4713 and 
He 11 4542 lines of the disentangled and observed spectra at the 
m aximum of separation and we compare these values with those 
of lConti & Alschuleii (119711) for stars of the same spectral types 
as both components. Moreover, we also measure the EWs of the 
above-quoted lines on synthetic spectra with the same spectral 
classifications as explained in Section |2l Even though the lu- 
minosity class is unknown for both components of HD 193443, 
the lines that we selected for this computation are only slightly 
affected by the luminosity class of the stars. Therefore, we de- 
cide to take into account the "canonical" EWs of all the stars 
with only similar subtypes. The ratios of these values provide 
a biightness ratio of about 3.9 + 0.4 between the piimary and 
the secondary components. The disentangled spectra, corrected 
from the brightness ratio, are shown in Fig. |3] 

The properties of both components thus show two stars with 
relatively close spectral types but with large mass and biightness 
ratios. This indicates that the primary is rather evolved relative 
to the secondary. First, we assume that HD 193443 is either an 
09III-H09.5V or an 09IH-09.5III system. In the former case, 
the com parison between the minimum masses and the values 
quoted in lMartins et al.l (I2005L 22 Mq and 16 for the primary 
and the secondary, respectively) suggests an inclination between 
17° and 21°. In the latter case, we obtain, from masses estimated 
to 30 Mq for the primary and 21 Mq for the secondary, an incli- 




-150 ^ ' < ' ' ' ' < ' ' ' ' < ' J 

0,5 1 

I'hase (0) 

Fig. 2. Top: He 1 447 1 , He ii 4542, He ii 4686, He 1 47 1 3 line pro- 
files of HD 193443 at different epochs. Bottom: RV curves of 
HD 193443 computed for a period of 7.467 days. Filled circles 
represent the primary RVs whilst the open circles correspond to 
the secondary. 



nation between 15° and 19°. Therefore, we consider that the in- 
clination of the system is very low and probably ranges between 
15° and 22°. From the formula of Eggleton ( 1983), we compute 
the projected Roche lobe radii (RRLsini) of 8.3 Rq and 5.6 Rq 
for the primary and the secondary components. These values in- 
dicate Roche lobe radii between 22 Rq and 32 Rq for the piimary 
and between 15 Rq and 22 Rq for the secondary. In both cases, 
the secondary (which would have a radius of 7 or 13 Rq de- 
pending on the case) and the primary (with a radius between 
13 Rq and 22 Rq) would not fill their Roche lobe. Secondly, if 
the stellar classifications were 09I for the primary and 09.5V 
for the secondary, the values of the inclination would be esti- 
mated between 15° and 21°, thereby implying Roche lobe radii 
between 23 Rq and 32 Rq for the primary and 16 Rq and 22 Rq 
for the secondary. In this case, the primary could fill its Roche 
lobe. At this stage and for the other SB2 systems, it is impor- 
tant to note that the analysis of the orbital parameters and the 
estimation of the inclination of the SB2 systems are only made 
to give a first idea of the geometry of the systems. However, 
we must be careful on the information provided here. Indeed, 
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Fig. 3. Disentangled spectra of the two components of 
HD 193443. The spectra are normalized accounting for the opti- 
cal brightness ratio of 3.9. The secondary spectrum is vertically 
shifted for clarity. 



the masses given in the tables of iMartins et al.l (l2005h are stan- 
dard values that also depend on other parameters such as e.g., 
logg, TefF and/or log(L/Lo). The values quoted in these tables 
can thus be different from the real parameters of both compo- 
nents. Therefore, it is common that the mass ratios computed 
from these standard values do not correspond to the real mass 
ratios determined from the orbital solutions. 

3.1.2. HD 228989 



Located in the Berkeley 86 young open cluster (lForbesl[l981b . 
this star was quot ed as an 09.5+09.5 binary system by 
iMassev et al.l (Il995h . However, no orbital solution of this object 
was yet determined. 

We collected eighteen observations of this object between 
Aug 2009 (HJD = 5049.9519) and Sep 2011 (HJD = 
5815.6460). These data clearly reveal, as it is shown in the up- 
per panel of Fig. |4] a binary system composed of two stars with 
rather similar spectral classifications. Our analysis of the refined 
RVs gives an orbital period of 1.77352 days. We compute from 
this value the orbital parameters reported in Table |3] We have 
computed two orbital solutions: an eccentric and a circular. The 
rms being smaller for the circular orbit, we decide to fix the 
eccentricity to zero. The RV curves computed with LOSP are 
shown in the lower panel of Fig. |4] As suggested from the disen- 
tangled spectra and the mass ratio, the two components seem to 
be quite close to each other in terms of the physical properties. 
From the raw (i.e., not corrected by the brightness ratio) disen- 
tangled spectra and the observed spectra, we determine spectral 
types of O 8.5 (logW = 0.238) and O 9.7 (log W = 0.791) for 
the primary and the secondary components, respectively. In ad- 
dition, we obtain, on the most deblended spectra, for the primary 
log W" = 0.027 and for the secondary log W" - 0.104, sug- 
gesting a main-sequence luminosity class for both components. 
Knowing the spectral classifications of these two stars, we es- 
timate a brightness ratio of about 1.2 + 0.1, the primary being 
slightly brighter than its companion. This brightness ratio cou- 
pled with the luminosities listed in tables of Martins et al. (2005) 
seem to agree with the luminosity classes of both components. 
Finally, we display in Fig. |5] the disentangled spectra corrected 
for the relative brightness. 

The comparison between the minimum masses and the 
masses published by Martins et al. (2005, 19 Mq and < 17 Mq 
for the primary and the secondary, respectively) suggests that 
the system would have an inclination between 42° and 50°. We 
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Fig. 4, Top: He 1 447 1 , He ii 4542, He ii 4686, He 1 47 1 3 line pro- 
files at different epochs. Bottom: RV curves of HD 228989 com- 
puted from Porb = 1.77352 days. Filled circles represent the pri- 
mary whilst the open circles correspond to the secondary. 



st ress however that the 09.7 subtype is not included in the table 
of lMartins et al.l (l2005l) . Therefore, we consider the values of the 
09.5 subtype as an upper limit. On the basis of Eggleton's for- 
mula (1983), we estimate RRLsini of about 5.7 Rq and 5.3 Rq 
for the primary and the secondary components. From the incli- 
nation range estimated from the masses, we compute a Roche 
lobe radius between 7.4 Rq and 8.4 Rq for the primary and be- 
tween 7.0 Rq and 8.0 Rq for the secondary. Such values seem to 
indicate that both components are very close to fill their Roche 
lobe (standard radii are of 7.9 Rq and < 1.2 Rq for the primary 
and the secondary, respectively). 

3.1.3. HD 229234 



Quoted as belonging to NGC6913 (iHumphrevsl Il978h. 
HD 22 9234 was classified as an 09If star bv ^ Massev et al.l 
([T99I w hilst iNegueruelal (l2004h reported an 09II classifica- 
tion and 'Morgan et al.l (11955) reported this sta r as 09.5III. 
Previous investigations of fLiu et al. (1989) and iBoeche et al.l 
(2004) proposed this star to be a binary system with an or- 
bital period of 3.5105 days. In addition, iBoeche et aH (|2004 
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4500 4600 4700 

x[A] 



Fig. 5. Disentangled spectra of the two components of 
HD 228989. The spectra are normalized accounting for a bright- 
ness ratio of 1 .2. The secondary spectrum is vertically shifted for 
clarity. 



published the first SBl orbital solution for this star. This SBl 
binary status was moreover confirmed by a recent analysis of 
iMalchenko & Tarasovl(l2009l) . 




4470 4540 4680 4710 

A[A] 



Fig. 6. Short-term variations of the Hei 4471, Hen 4542, 
Hen 4686 and He: 4713 fine profiles of HD 229234. The rest 
wavelengths are given by the dashed lines. 



We obtained nineteen spectra of HD 229234 between Sep 
2008 (HJD = 4717.3892) and Aug 2010 (HJD = 5421.4347). 
This object obviously displays significant variations in its RVs 
on a short timescale. Fig. |6] shows the Doppler shifts of several 
spectral lines for six consecutive nights. The periodic motion is 
clearly visible for all the spectral lines. The mean radial veloci- 
ties (RVs) ± the standard deviation computed from the RVs are 
-8.7 + 33.8, 4.9 + 34.4, 1.7 + 36.5 and -6.5 + 34.0 km s"' for 
the He 1 447 1 , He ii 4542, He n 4686 and He 1 47 1 3 lines, respec- 
tively. No SB2 signature was found in the optical data. We de- 
termine on the basis of the HMM technique an orbital period of 
3. 510595 + 0.000245 d ays, agreeing with th e previous estimates 
of iLiu et al.1 11989) and lBoeche etal.1 ( |2004 . The parameters of 
the SBl orbital solution are listed in Tableland the RV curve is 



displayed in Fig.[71 These parameters are globally close to those 
obtained bv Boech e et al.l (l2004l) . With the absence of uncertain- 
ties in the orbital solution given by these authors, the comparison 
remains however diflicult. The main diff'erence is regardless ob- 
served in To value because of a different definition of the phase 
= 0.0. 

We also derive, for the main component of this SB 1 system, 
an O 9III classification from the different EW ratios log W = 
0.372, logW" = 0.280 and logW" = 5.325. This spectral 
type appears to be in agreemen t with the different cl assifications 
found in the literature. Unlike iMassev et al.l (1 19951) . our classi- 
fication is not only based on a visual comparison with spectra 
of the atlas of IWalbom & Fitzpatrickl (11990) . but also relies on 

quantitative values. 

From the mass of the primary (22.0 + 3.0 Mm iMartins et al.l 

l2005h and the /mass value given in Table [3] we determine, as 
minimum mass for the secondary component of the system, a 
value of about 3.0 + 0.3 Mq. 



4 - 




0.5 1 

Fhase (0) 

Fig. 7. RV curve of HD 229234 computed from Poa = 
3.510595 days. 



3.2. Presumably single stars 
3.2.1. HD 193514 

This star was classified as 07Ib(f) bv ' Repolust et all ( |2004 . 
These authors also reported on a weak emission in the He n 4686 
line profile. Chan ges from day to day were already observed by 
lUnderhill (Il995 h in the Ha profile. However, none of these au- 
thors did attribute these variations to a binary nature. 

The twenty-two spectra of HD 1935 14 taken from Aug 2004 
(HJD = 3246.5202) to Aug 2010 (HJD = 5420.4585) do not 
confirm any significant RV variations. The upper parts of Fig. [8] 
show that the Doppler shifts of the spectral lines are rather small 
(undetectable by eye) on several (short and long) timescales. 
We indeed determine RVs = -18.0 + 5.2,-7.0 + 5.0 and 
-15.6 ± 5.0 km S-' for the He i 4471, Hen 4542 and He i 4713 
lines. In the lower parts of Fig. |8] we show the TVS computed 
for these different spectral lines. The variations above the dot- 
ted lines are considered as significant under a significance level 
of 1%. Therefore, we see that the TVS of the Hen 4542 and 
the He I 4713 line profiles do not show any substantial varia- 
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Fig. 8. The He i 4471, Hen 4542, Hen 4686 and He i 4713 line 
profiles of HD 193514 according to HJD. The rest wavelengths 
are given by the dashed lines. At the bottom, TVS of each spec- 
tral line computed from the Aurelie data. The dotted line illus- 
trates the 1 % significance level for the variab ility evaluated fol- 
lowing the approach of .FuUerton et al] (Il996l) . 



tions, though the Hei 4471 and especially the Hen 4686 lines 
do. However, since this variability is not detected on all the line 
profiles, we conclude that it is unlikely to be due to binarity but 
rather arises in the stellar wind itself. However, we do not have 
enough echelle s pectra to characterize the changes reported by 
lUnderhilH (Il995h in the line profile of Ha. We therefore consider 
HD 193514 as a presumably single star. 

From the echelle spectra, we determine log W' - 0.014 and 
\ogW = 0.253, yielding for HD 193514 a spectral type 07- 
07.5 and a luminosity class similar to a giant (111). The strong 
Niii 4634-41 lines coupled with the Hen 4686 line (for which 
the emission component is similar to the absorption one) sug- 
gest an (f) suffix, which agrees with a giant luminosity class. 
Therefore, we classify HD 193514 as an 07-07.5 111(f) star 

3.2.2. HD 193595 

HD 193595 is locat ed in the Berkeley 86 y oung open cluster 
(iMassev etalJll995h . Classified as an O 6V (iGarmanv & Vacca 
[l99ir this object was poorly investigated in terms of RV. 

We measure the RVs on nine spectra taken at different 
epochs, between Sep 2008 (HJD = 4711.3816) and Aug 2010 
(HJD = 5421.4977). These data reveal rather constant RVs 
for this star: RVs = -13.7 + 2.2, -8.1 + 2.8, -4.5 + 2.6 and 
-12.5 + 1.7 km s-i for the He 1 4471, Hen 4542, Hen 4686 and 
He I 4713 lines, respectively. We detect, as shown in the lower 
panels of Fig.|9] low amplitude variations (at a significance level 
of 1%) only in the TVS spectrum of the He i 4471 line. The re- 
sults of the TVS associated to the RV measurements thus seem 
to imply that the observed variations are probably not due to bi- 
narity. By assuming a single status for this object, we derive for 
HD 193595 an 07V steUar classification (log W = -0.036 and 
log W = 0.052). 



471 1,3816 
4712. 4309 




Fig. 9. Same as Fig.[8]but forHD 193595. 



3.2.3. HD 193682 

This star is the hottest object in our sample. HD 193682 seems 
to have a spectral type oscillating between 05 (lHiltneiiri956h 
and 04 (04111(f), Gai-manv & Vacca 1991). We obtained ten 
spectra of this object between Sep 2008 (HJD = 4711.4038) 
and Aug 2010 (HJD = 5418.5280). These data show rather 
broad line profiles, indicating a rather large projected rotational 
velocity. We determine from the Fourier transform method of 
Simon-Diaz & Herrero (2007) a v sin / = 150 km s"' for this ob- 
ject. This moderately fast rotation yields broader lines that pro- 
duce larger uncertainties on their centroid and hence poorer ac- 
curacy of the RVs. However, we do not perform a least-square 
fit with a synthetic profile broadened by rotation to determine 
the RVs of this star because the shapes of the strongest lines 
are still Gaussian. We determine RVs of about -48.6 + 11.3, 
-39.7 + 7.7 and -55 . 1 + 1 1 . 1 km s- ' for the He 1 447 1 , He n 4542 
and Hen 4686 lines, respectively. Such a high 1-cr dispersion is 
probably due to the line broadening rather than to binary motion 
because this RV dispersion is smaller on the strongest spectral 
lines. Furthermore, we do not detect any secondary signature 
in the observed spectra nor any change of the line width as a 
function of time, as it is shown in the upper panels of Fig. [TO] 
Indeed, we see from this figure that the spectral lines have al- 
ways their centre at the same side of the dashed lines, which 
represent the rest wavelengths. We compute the TVS for these 
spectral lines (lower panels of Fig. [TOl i and we find a signifi- 
cant (although weak) variation in the line profile of He n 4686, a 
spectral line generally affected by strong stellar winds. We can- 
not consider these changes as due to binarity because they are 
not detected in all the spectral lines, notably in the Hei 4471 
and the He n 4542 lines, which suggest that these variations are 
probably intrinsic to the stellar wind rather than attributable to a 
putative companion. 

We then calculate, from the Hei 4471-Hen 4542 ratio, 
log W' = -0.5 14, corresponding to an O 5 subtype. Furthermore, 
Walborn's criterion suggests to add a (f) suffix, indicating a mod- 
erate emission for the N ni 4634-4 1 lines and a moderate absorp- 
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Fig. 1 1. The He 1 447 1 , He ii 4542, He ii 4686 and He 1 47 1 3 line 
profiles of HD 194094. 



tion for the Hen 4686 line, similar to giant (111) stars. To sum- 
marize, we classify HD 193682 as an 05111(f) star. This spectral 
classification is similar to those found in the literature. 



3.2.4. HD 194094 

During our spectroscopic campaign, HD 194094 was only ob- 
served 5 times: twice in Jun 2009, twice in Aug 2010 and once 
in Sep 201 1 . These data however allow us to investigate RV vari- 
ations on timescales from one week to one year. These five spec- 
tra do not exhibit any evidence of a putative companion on such 
timescales nor any periodic motion relative to the rest wave- 
lengths (dashed lines in Fig. [TTI ). We measure the RVs on the 
He I 4471, Hen 4542, Hen 4686 and He i 4713 Hnes, reporting 
Ws = -18.2+3.9,-11.3+3.8, -9.7+8.9 and -13.9+6.1 kms"', 
respectively. The standard deviations on these measurements in- 
dicate that HD 194094 is probably single and that the RV vari- 
ations detected on the He ii 4686 line are likely due to a vari- 
able stellar wind. However, unlike the stars previously analysed, 
the Aurelie data for HD 194094 are not sufficiently numerous 
nor sufficiently spread over the entire campaign to compute the 
TVS. We also cannot include the SPM data in the computa- 
tion of the TVS, as we would need to convolve the data so that 
they all have the same resolution. Indeed, unlike HD 46150 in 
iMahv et al.'(2009'). we do not have, in the [4450-4900] A wave- 
length domain, diagnostic lines similar to the interstellar Nai 
lines at 5890 A and 5896 A to compare the convolved spectra. 

Finally, by assuming that this object is single, we com- 
pute, from an echelle spectrum of HD 194094, log W - 0.232, 
log W" = 0.174 and log W" = 5.204, yielding an O 8.5111 star 
This spectral type is slightly earlier than reported in the literature 
(O 9111, ^Hiltner 1956). However, the error-bars on the determi- 
nation of the spectral type do not allow us to claim that the new 
classification is significantly different from the older ones. 



3.2.5. HD 194280 

Presented as a prototype of the OC9.71ab class, HD 19428 ap- 
pears to be a carbon-rich star wit h a depletion in nitrogen (iGovl 
I1973I: IWalbom & Howarthll2000h . According to these authors, 
the profiles of Hei 5876 and Ha would indicate that the stellar 
wind is weak. 
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Fig. 12. Sample of spectra taken from the entire dataset. This 
figure is similar to Fig. |8]but for HD 194280. 

We obtained forty-two spectra of this target with the Aurelie 
spectrograph in the [4450 - 4900] A wavelength domain be- 
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tween Oct 2004 (HJD = 3286.3646) and Nov 2007 (HJD = 
4423.2513). The spectral lines sampled over different timescales 
are displayed in the upper panels of Fig. [12] This figure shows 
that the Doppler shifts of the lines are small, which is confirmed 
by the RVs. These values for the He 1 447 1 , He ii 4542, He ii 4686 
and He 1 47 1 3 lines are indeed equal to -24.9 ± 4.2, - 1 1 .6 + 7 .6, 
-5.4 + 5.0 and -17.1 ± 4.8 km s"', respectively, suggesting that 
this star is presumably single. From the He 1 4471-Heii 4542 ra- 
tio, we obtain log W - 0.742, which agrees with the 09.7 sub- 
type. Furthermore, the object presents strong line-profile varia- 
tions notably in the Hei 4471 line. The results of the TVS are 
shown in the lower panels of Fig. [12] Since the profile of the 
He n 4686 line does not appear in emission and since the profiles 
of the He I 5 876 and Ha Unes did not seem to indicate strong 
stellar wind (lGovlll973 ^. these variations are unlikely due to 
the wind. Therefore, another possible explanation for this strong 
variability in He i 4471 could be pulsations. 

3.2.6. HD 228841 



iMassev et al.l ( Il995h located this object in the Berkeley 86 young 
open cluster, along with HD 193595 and HD 228989. This 
sta r was success ively classified as 07y ((f)) and 06.5Vn((f)) 
by iMassev et al.l (1995) and So ta et all f2011, and references 
therein). 

Eight optical spectra were collected from Sep 2008 (HJD - 
4712.3982) to Aug 2010 (HJD = 5418.5980), mainly in the 
[4450 - 4490] A range. The spectra show broad lines which 
means that HD 228841 is likely a rapid rotator We indeed 
see, in the upper panels of Fig. [13] that the spectral line 
widths remain constant as a function of time, which does 
not seem to support the presence of a secondary component. 
Therefore, we determine fro m the Fourier transform method of 
ISimon-Diaz & Herrerol (I2007J the projected rotational velocity 
of the star. We obtain a value of 3 17 km s ' . The mean RVs de- 
termined by least-square fit between a rotation profile and the 
spectral lines amount to RVs - -48.4 + 8.9, -25.5 ± 6.5 and 
-44.9 ± 25.3 km s-i for the He 1 4471, He ii 4542 and He 1 4713 
lines, respectively. Given the v sin /, the standard deviations gen- 
erally smaller than 15 km s"' favour a single status for this ob- 
ject. In addition, the variations in the line profiles computed with 
the TVS (see the lower panels of Fig. \T3[ are not significant, 
which agrees with the presumably single status for HD 228841. 
The high dispersion detected for the He i 4713 line is probably 
due to the weakness of this line. 

We then derive the spectral type of this star from Conti's cri- 
terion, log W is equal to -0.01 1, wh i ch ind icates an O 7 star, as 
already determined by IMassev et al.l (Il995l) . We stress however 
that log W" could not be estimated given the broadness of the 
lines and thus the severe blend of the Si iv 4089 line with H6. 

4. 0-type stars in the Cyg 0B3 association 

4. 1 . Presumably single stars 
4.1.1. HD 190864 



Quoted as belonging to NGC6871 (iHumphrevsl 1 19781) . 
HD 190864 is one of the brightest O-type stars of CygOB3. The 
spectral class ifications d erived for this objec t range between O 6 
(lHiltner'1956) and O 7 (IConti & Leeplll97 4) but agree on a gi- 
ant luminosity class. Previous investigations (e.g., Plaskett 1924) 
quoted this star as SBl but no solution for its orbital motion 
has been published. iMeisell d 19681) listed this star as visual dou- 
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Fig. 13. Same as Fig.lbut for HD 228841. 



ble star with Am - 0.5 (O6+B0.5Vp). Although these authors 
did not mention the angular separation, an B 0.5Vp spectral type 
could match with the spectral ty pe of HD 227586 located a t 1' 
from HD 190864. More recendv. lGarmanv & St encei'('1992') at- 
tributed to HD 190864 a spectral type of B III-hO 6.5111(f), as- 
signing it an SB2 status. However, it is possibl e that the B II I 
component is the visual component reported bv lMeisel (1 19681) . 
Indeed, the latter author also mentioned that the estimate of the 
luminosity class of the B component was uncertain and could 
be higher, thereby cou l d be attributed to the B 111 component 
of lGarmanv & Stencell (Il992h. From seven observations spread 
over two years. Underbill did however not detect any 

significant RV variations. According to this author, it is un- 
likely that the differen ce between the projected rotational ve- 
locity given by j Conti & EbbetsI (Il977i, 69 km s ') and that of 
iHerrero et all (Il992i 105 km s ') is due to a possible varying 
blend of the profiles of two putative components but rather seems 
to originate from an emission component in some absorption 
lines. 

We obtained seventeen spectra between Aug 2001 (HJD = 
2134.4733) and Aug 2010 (HJD = 5419.3479). This dataset 
does not show any significant RV changes, except for the 
He II 4686 line. The mean RV and the standard deviation mea- 
sured on this line are indeed of RV - 4.6 + 11.5 km s ' 
whilst the other values deri ved on the He i 4471, He ii 4542 and 
Hei 4713 lines amount to RVs = -15.5 + 2.5, -6.3 + 2.4 and 
-16.5 + 3.8 km s"', respectively. Since our observations cover 
both short and long timescales, we can exclude that HD 190864 
is a binary system with an orbital period shorter than about 3000 
days. From an intense survey devoted to HD 190864, we have 
never detected the signature of a putative companion. If a B com- 
ponent was pres ent in the spectrum of H D 190864, we should see 
it. According to lSchmidt-Kaleil(ll982l) . a B 0.5V or a B III com- 
ponent should be only 6 times fainter than an O 6.5111(f) star and 
the mass ratio should be between 1 and 4. Therefore, such a com- 
ponent should be visible in the observed spectra, which is not the 
case. It is thus unli kely that the differe nce b etween the two value s 
of V sin /, given by IConti et all (Il977h and iHerrero et all (Il992l) . 
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Fig. 14. Same as Fig.|8]but for HD 190864. 



originates from binarity. Fuithermore, if we co nsider that the RV 
published by Conti et al. (-15.3 km s~ Ml977b was measured on 
a He I line, this could constitute an additional clue that the star 
is presumably single. This thus reveals the importance of a ho- 
mogeneous dataset to study the binary fraction of massive star 
populations (as mentioned in Sect.|2]i. 




4840 4850 4860 4870 4880 
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Fig. 15. Mean spectrum and TVS of HD 190864 computed from 
the Aurelie data of the H/3 line. The dotted line illustrates the 
1 % significa nce level for the vari ability evaluated following the 
approach of iFullerton et al.l(ll996l) . 

Though no signature of a secondary companion is found in 
the HD 190864 spectrum, the TVS analysis reveals significant 
variations of some line profiles, mainly in the red wings of the 
He 11 4686 and H/3 lines (lower panels of Figs. [14] and Fig.fTsTl. It 
is however unlikely that these variations are due to the existence 



of a secondary star. If that was the case, the variability pattern 
should be once again observed in all the spectral lines which is 
clearly not the case. Therefore, we classify this object as pre- 
sumably single. 

Finally, we derive from Conti's criterion a spectral classifi- 
cation of 06.5 (log W = -0.193), consistent with the previous 
estimates. The presence of moderate emission in the N iii 4634- 
41 lines and a rather weak Hen 4686 line suggest to add the 
(f) suflix. Therefore, we derive that the spectral classification for 
HD 190864 is 06.5111(f). 

4.1.2. HD227018 

Located at about 1.2° south-west of HD 190864, HD 227018 is 
rather faint {V - 8.98) for our survey. While its spectral type is 
relatively well known (06.5-07), its luminosity class is poorly 
constrained. Indeed, Herrero et al. (1992) quoted this star as an 
O 6.5111 while iMassev et aL (,1995.) attributed it an 7V((f)) 
type. 
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Fig. 16. Same as Fig.[8]but for HD 227018. 



Between Sep and Oct 2008, we collected eight spectra with 
the Aurelie spectrograph and we increased our dataset with 
two additional spectra taken in August 2010 with the same in- 
strument. Th e me asurements of the RVs reveal rather constant 
values with RVs = 15.6 + 3.4, 21.7 + 4.0, 23.8 + 3.3 and 
16.6 + 3.9 km s ' for the Hei 4471, Hen 4542, Hen 4686 and 
He I 4713 lines, respectively. We see from the upper panels of 
Fig.[T6]that the Doppler shifts of the spectral lines are relatively 
small relative to the rest wavelengths. Moreover, the TVS spec- 
tra (lower panels of Fig. [T6t only exhibit small variations in the 
He I 447 1 line at the limit of being significant while the other 
lines are relatively stable, thereby indicating that HD 227018 is 
presumably single. In order to derive the spectral type of this 
star, we compute, from the Hei-Heii ratio, log W = -0.201 
which corresponds to an 06.5 star We add the ((f)) suffix rem- 
iniscent of a main-sequence star because the object shows weak 
Niii 4634^1 emissions and strong Hen 4686 absorption lines. 
Therefore, we classify this star as 06.5V((f)). 
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4.1.3. HD 227245 

HD 227245 is known as an O 7V jGarmanv & Vaccalll99ll) but 
its multiplicity has never been investigated until now. From five 
observations taken between Aug 2009 (HJD - 5049.8148) and 
Jun 2011 (HJD = 5727.7497), we determine RVs = 5.9 + 5.4, 
9.8 + 3.9, 10.8 + 6.2 and 6.0 + 7.3 for the He 1 4471, Hen 4542, 
Hen 4686 and Hei 4713 lines, respectively. Figure [TT] shows 
these spectral lines at different HJD. We stress that a normaliza- 
tion problem occurs at HJD = 5723.9405 for the He i 4471 line 
but it does not affect the determination of the RV for this line. 
The standard deviations computed from the RVs are not consid- 
ered significant because they are generally smaller (except for 
He I 4713) than the 7-8 km s ' threshold defined as our vari- 
ability criterion. The exception for the He 1 4713 line can be ex- 
plained by the weakness of this line and the signal-to-noise ratio 
of the data. HD 227245 is thus considered as presumably single. 
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Fig. 17. The He 1 447 1 , He ii 4542, He ii 4686 and He 1 47 1 3 line 
profiles of HD 227245. 



(see Fig. [HI). Indeed, we measure RVs = -26.8+4.6, -26.2+4.5, 
-28.7 + 2.9 and -30.0 + 4.3 km s"' for the He 1 4471, He ii 4542, 
Hen 4686 and Hei 4713 lines, respectively. We thus consider 
the star as presumably single. 
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Fig. 18. The He 1 447 1 , He ii 4542, He ii 4686 and He 1 47 1 3 line 
profiles ofHD 227757. 



The spectrum of HD 227757 displays a stronger He i 4471 
line than the Hen 4-542 line, indicating a spectral type 
(IWalborn & Fitzpatrickl 19901) later than O 7. Due to the late sub- 
type and the small v sin / of the star (~ 45 km s ' ), we can clearly 
distinguish an On line blended with the Hei 4471 line. By re- 
moving the contribution of the former line (i.e., by fitting two 
Gaussian profiles, one on the O ii line and one on the He i line), 
we measure log W - 0.332, corresponding to an O 9 star. We 
also obtain log W" = 0.081 and log W" = 5.485 which indi- 
cate a main-sequence luminosity class (V). Therefore, we clas- 
sify HD 227757 as09V. 



The spectrum of HD 227245 exhibits He ii slightly weaker 
than He i, indicating a mid 0-type star. We obtain log W - 
-0.034 and log W = 0.218, corresponding to a spectral type 
07III. However, the weak Niii 4634^1 lines coupled with a 
strong He ii 4686 line suggest an ((f)) suffix, indicating rather a 
main-sequence luminosity class. This difference could come ei- 
ther from the definition of the continuum level in the normalized 
spectra or because Conti's criterion is barely applicable to this 
spectral type. Therefore, we classify HD 227245 as an O 7V- 
III((f)) s tar. This spectral classifica tion is in agreement with that 
given bv lGarmanv & Vaccal (1199 lb . 

4.1.4. HD 227757 

HD 22 775 7 (V ^ 9.25) was r eported by iGarmanv & Stencell 
(Il992h and lHerrero et all (11992*) to be an 09.5V star. From six 
spectra taken between Jun 2009 (HJD = 4988.8797) and Jun 
2011 (HJD - 5728.7408), we do not observe significant varia- 
tions in the RVs of HD 227757 on short nor on long timescales 



5. 0-type stars in the Cyg 0B8 association 

5.1. Presumably single stars 
5.1.1. HD191423 

Better known under the name of Howarth's star, HD 191423 
is the fastest rotator known to date among Galactic O stars 
and it was quoted as the prototype of the ONn stars (nitrogen- 
rich O star with broad diffuse lines, ,Walborn ,200 3). Its pro- 
jected rotati£nalvek)citY_was estimated to be around 400 km s ' 
(e.g., iPennv & GiesI l2009l computed a vsin/ in the range of 
[336-436]kms-i). 

We obtained nine spectra of this star between Aug 2004 
(HJD = 3247.3339) and Aug 2 010 (HJD = 5421.5876 ) From 
the Fourier transform method ofi imon- Diaz & Herrerol(l2007h . 
we determine v sin / = 410 km s"', in good agreement with pre- 
vious estimates. Moreover, we see in the time series of the spec- 
tra (shown in upper panels of Fig. [19]) that the line widths remain 
constant as a function of time, thereby suggesting rather a rapid 
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Fig. 19. Same as Fig.|8]but for HD 191423. 
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Fig. 20. Same as Fig.lbut for HD 191978. 



rotator than a binary system. The RVs measured by least-square 
fit give us average values of RVs - -44.0 + 3.0, -58.8 + 18.4, 
-48 .7 + 5 .0 and -56.5 + 4.8 km s- ■ for the He 1 447 1 , He n 4542, 
Hen 4686 and He i 4713 lines, respectively. The RV uncertain- 
ties for rapid rotators are typically larger because the centroids 
of the lines are not well defined. This could explain the large 
standard deviation observed for the weak He n 4542 line whilst 
the other lines are rather stable in RV, thereby indicating a pre- 
sumably single star Moreover, after having computed the TVS 
spectra for these characteristic lines only from the Aurelie data 
(lower panels of Fig. \T% . we detect no significant variation 
of line profiles, which agrees with the likely single nature of 
HD 191423. 

Therefore, we compute log W' = 0.396 and log W" - 
5.187, assigning to this star an 09111 type. We are not able to 
estimate log W" given the blend of the Si iv 4089 line with H6. 
We also add the suffixes n and N, because of the broadness of 
the lines as well as the strong relative intensity of the nitro- 
gen lines, thereby giving an ON 9111n spectral classification for 
HD 191423. 

5.1.2. HD 191978 



HD 191978 was classified as an O 8 star bvlGovl (Il973l) . The 
investigation of its RVs bv lAbt&Big"gsl(ll972r did not reveal 
any significant variations. 

We collected eleven spectra of HD 191978 over a timescale 
of about 2170 days between Aug 2004 (HJD = 3247.4196) 
and Aug 2010 (HJD = 5419.4729). The measurements made 
on these data show relatively constant RVs. We can see in the 
upper panels of Fig. |20]that no clear Doppler shift exists relative 
to the rest wavelengths of four spectral lines. The mean RVs and 
their standard deviations are RV - -18.1 + 3.9 km s ' for the 
He I 4471 line, RV = -10.4 + 5.3 km s ' for the Heii 4542 
line, RV - -9.5 + 5.0 km s ' for the Hen 4686 line, and 
RV = -11.8 ± 6.2 km s ' for the Hei 4713 line. However, 
the TVS spectra (lower panels of Fig. l20t indicate significant 
line profile variations for the Hei 4471 and marginal ones for 



the He ii 4686 line, whilst insignificant line profile variations are 
found for the Hen 4542 and Hei 4713 lines. Since these vari- 
ations are not detected in all the line profiles, it is not likely 
that they are linked to binarity. Therefore, we consider this star 
as presumably single. The determinations of the EWs for the 
diagnostic lines yield log W = 0.120, logW" = 0.227 and 
log W" = 5.398, which correspond to an O 8111 star. 

5.1.3. HD193117 

As several O stars in our sample, HD 193117 has received lit- 
tle attention in the past years. We obtained eight spectra be- 
tween Sep 2008 (HJD = 4711.524 5) and Aug 2010 (HJD = 
5420.6211). Quoted as an O 9.511 by Govl ([T973ft . this star does 
not present any significant variations in its RVs. This can be 
seen in the upper panel of Fig. [2T| We indeed measure RVs = 
-29.1 + 2.4, -16.6 + 7.1, -13.5 + 4.8 and -26.0 + 4.7 km s"' 
for the He i 4471, He ii 4542, He ii 4686 and He 1 4713 Unes, re- 
spectively. In addition to these results, the TVS (lower panels of 
Fig. |2TJ shows no significant change in the line profiles of the 
star. Since we find no direct evidence of variations in the spec- 
trum of HD 1931 17, we report this star as presumably single. 

The diagnostic line ratios, Hei 4471-Heii 4542 and 
Si IV 4089-Hei 4143, give log W = 0.343 and log W" = 0.256. 
These results correspond to an O 9111 star. Furthermore, we com- 
puted log W" = 5.155 which agrees with this spectral classifi- 
cation. 



6. 0-type stars in the Cyg 0B9 association 

6.1. Gravitationally bound systems 
6.1.1. HD 194649 

Often q uoted as an 06.5 star, HD 194649 was reported by 
'Mulled (11954 to be a binary system. However, no orbital so- 
lution is found in the literature. Our dataset is composed of sev- 
enteen spectra taken between Jun 2009 (HJD = 4992.8930) and 
Sep 2011 (HJD = 5820.6424). These data reveal a clear SB2 
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Fig. 21. Same as Fig. |8]but for HD 193 1 17. 



signature with a secondary component moving in anti-phase rel- 
ative to the primary. This motion is clearly observed over two 
consecutive nights, indicating a short-period binary system. By 
applying the Fourier method of HMM to the RVs refined by 
the disentangling programme, we derive an orbital period of 
3.39294 + 0.00139 days. As we did for the previous binary sys- 
tems, we compute two orbital solutions: an eccentric one and 
a circular one. However, the best orbital solution (i.e., with the 
smallest rms) is achieved for a circular orbit. The orbital param- 
eters from this solution are given in Tableland the RV curves 
are shown in Fig. |22] The primary appears to be significantly 
more massive than the secondary. Moreover, we have a differ- 
ence of about 16 km s"' between the systemic velocities of both 
components. This discrepancy could indicate that the secondary 
lines are more affected by the velocity field of a strong stellar 
wind, although such an interpretation is at odds with the spectral 
types derived below which suggest an earlier and more evolved 
primary star. 

We measure the EWs to determine the spectral classifica- 
tion on the disentangled and observed spectra. We find log W - 
-0.267 and log W - 0.124 for the primary and secondary com- 
ponents, respectively, corresponding to 06 and 8 subtypes. 
For the primary, it is not possible to compute log W" because 
the star is too early for Conti's criterion. Therefore, by focusing 
on the disentangled spectra (not yet corrected for the brightness 
ratio), the rather moderate Nin 4634-41 and Hen 4686 lines 
suggest to add an (f) tag to the primary spectral type, thus sug- 
gesting a giant luminosity class. For the secondary, we determine 
log W" - -0.03, corresponding to a main-sequence star (V). We 
thus conclude that the spectral classifications for both stars are 
06111(f) and 8V for the primary and the secondary, respec- 
tively. On the basis of these spectral classifications, we derive 
a brightness ratio of about 4.7 + 0.4, which agrees with a pri- 
mary star more evolved than its companion. Finally, we use this 
brightness ratio to correct the individual spectra obtained by dis- 
entangling. The resulting spectra are displayed in Fig.l23l 

From such spectral classifications, the tables of lMartins et al.l 

(l2005i) indicate radii and masses of 15 Rq and 35 Mq for an 




0,5 



I'hase (0) 

Fig. 22. Top: Hei 4471, Hen 4542, Heii 4686, He i 4713 line 
profiles at different epochs. Bottom: RV curves of HD 194649 
computed from f orb = 3.39294 days. Filled circles represent the 
primary whilst the open circles correspond to the secondary. 



06111(f) star and of 8 Re and 21 Mq for an 8V star From 
these masses and the minimum values provided in Table |3] we 
compute an inclination between 26° and 32 ° for the sy stern . 
Furthermore, we determine, on the basis of lEggleton] (Il983l) . 
RRL sin i of 8.3 Rq and 5.4 Rq for the primary and the sec- 
ondary, respectively, thereby giving a radius between 16 Rq and 
19 for the primary Roche lobe and between 10 Rq and 12 Rq 
for the secondary one. If the standard values are close to those 
of both stars, it thus appears unlikely that both components of 
HD 194649 fill their Roche lobe. 

6.2. Presumably single stars 
6.2.1. HD 194334 

HD 194334 was classified as an O 7.5V star by |Go3 (Il973h . We 
collected eleven spectra between Sep 2008 (HJD = 471 1.5581) 
and Aug 2010 (HJD = 5420.4943). We compute Ws = 
-15.3 ± 4.6, -3.8 ± 5.2, 23.1 + 12.7 and -15.6 + 5.3 km s"', 
for the He i 447 1 , He ii 4542, He ii 4686 and He 1 47 1 3 Unes, re- 
spectively. These velocities reveal a significant variability in the 
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Fig. 23. Disentangled spectra of the two components of 
HD 194649. The spectra are normalized accounting for the 
brightness ratio of 4.7. The secondary spectrum is vertically 
shifted for clarity. 



6.2.2. HD 195213 

Reported as an O 7 star by iGovl (Il973l) . HD 195213 has been 
observed ten times during our campaign between Sep 2008 
(HJD = 4711.5962) and Aug 2010 (HJD = 5420.5309). Our 
data do not reveal any significant variations of the Doppler shifts 
for the main lines (see upper panels of Fig.l25ll. We indeed com- 
pute RVs = 2.1 + 4.5, 9.6 + 4.6 and -2.0 + 4.2 km s"' for 
the He I 4471, Hen 4542 and He i 4713 lines, respectively. The 
He II 4686 line exhibits a pattern with bf a strong emission in 
the line core that almost completely masks the absorption pro- 
file, thereby probably indicating that the stellar wind of the star 
is relatively strong. We detect, through the TVS analysis (lower 
panels of Fig. |25] ), significant variations for the Hei 4471 and 
Hen 4686 lines. Since the TVS does not indicate line profile 
variations in all the lines and since the RV variations are not sig- 
nificant, we assume that the observed variability is likely due to 
the stellar wind rather than to a putative companion. 



He II 4686 line while the other lines are rather constant (upper 
panels of Fig.l24ll. The mean RV of the He ii 4686 line is positive 
whilst the others are negative. This could suggest that this line 
is formed in the wind. Moreover, changes in certain line profiles 
are also detected in the TVS spectra (lower panels of Fig. l24t . 
as notably for the Hei 4471 and Hen 4686 lines. However, no 
similar change is observed in the He ii 4542 and He 1 4713 lines, 
which strengthens the assumption that these variations are likely 
produced in the stellar wind. By assuming that HD 194334 is a 
single star, we determine log W - -0.002 and log W" = 0.294, 
thus giving for HD 194334 an O 7-7.5111 spectral classifica- 
tion. The luminosity class is confirmed by Walborn's criterion. 
Indeed, the spectrum of HD 194334 shows moderate N iii 4634- 
41 emissions and He ii 4686 absorption lines, thereby suggesting 
a giant luminosity class. Therefore, HD 194334 is classified as 
an 07-0 7.5111(f) star. 




4470 4540 4680 4710 



A|A] 

Fig. 24. Same as Fig.|8]but for HD 194334. 




4470 4540 4680 4710 



A[A] 

Fig. 25. Same as Fig.[8]but for HD 195213. 



From the echelle spectrum, we measure log W - -0.026, 
corresponding to an O 7 star and log W" = 0.264, indicating a 
giant luminosity class (111). Moreover, we add the (f) suffix to 
the spectral type because the Niii 4634^1 emissions are strong 
whilst the He ii 4686 absorption is rather weak. We thus suggest 
an O 7111(f) spectral classification for this star 

7. Discussion 

7.1. Observational biases 

We performed Monte-Carlo simulations to estimate the proba- 
bility to detect binary systems on the basis of our temporal sam- 
pling, for the 15 presumably single stars. For that purpose, we 
randomly draw the orbital parameters of 100000 binary systems. 
We consider that these systems are not detected if the standard 
deviation of the RVs is smaller than 7-8 km s"' to be consis- 
tent with our binarity criterion. For the fast rotators, we request 
standard deviations of at least 15 km s"' . The period distribution 
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Table 3. Orbital solutions for the four newly-detected binary systems 



Parameters 


HD 193443 


HD 228989 


HD 229234 


HD 194649 




Primary Secondary 


Primary Secondary 


Primary 


Primary Secondary 


P[d] 


7.467 ± 0.003 


1.77352 ± 0.00041 


3.51059 ±0.00175 


3.39294 ± 0.00139 


e 


0.315 ±0.024 


0.0 


0.0 


0.0 




289.7 ± 5.8 








To (HJD) 


3246.120 ± 0.106 


5047.771 ± 0.003 


4714.326 ± 0.01 1 


4991.082 ± 0.009 


7 [km s"'] 


-10.0 ±1.2 -16.6 ±2.0 


-3.7 ±1.8 -7.8 ±1.9 


-16.3 ±0.8 


-13.8 ±2.1 -29.5 ±3.1 


a: [km s-'] 


38.8 ±1.5 91.3 ±3.4 


194.3 ± 2.5 220.3 ± 2.9 


48.5 ± 1.1 


76.3 ±2.5 192.1 ±6.2 


a sin i [/?©] 


5.4 ±0.2 12.8 ±0.5 


6.8 ±0.1 7.7 ±0.1 


3.4 ±0.1 


5.1 ±0.2 12.9 ±0.4 


M sin^ (■ [Mo] 


1.0 ±0.1 0.4 ±0.1 


7.0 ±0.2 6.1 ±0.2 




4.9 ±0.4 1.9 ±0.1 


q m^lM^) 


2.35 ±0.11 


1.13 ± 0.02 




2.52 ±0.13 


/mass [Mo] 






0.042 ± 0.003 




rms [km s"'] 


5.08 


8.37 


3.04 


11.13 



Notes. The l-cr error-bars on the parameters are provided by LOSP, except for the orbital period whose uncertainty is determined on the basis of 
the natural width of the peaks in the Fourier power spectrum. The orbital solution of HD 193443 was computed assuming that the uncertainties on 
the secondary RVs are three times as large as for the primary. 



is chosen to be bi-uniform in log scale as in lSana et al ] (l20Q9h . 
Under this assumption, log P is following a bi-uniform distribu- 
tion between 0.3 and 1 .0 for 60% of the systems and between 
1.0 and 3.5 for the remaining 40%. The eccentricities are se- 
lected uniformly between 0.0 and 0.9. As Rau w et al.l (1201 ll) al- 
ready did, we assume that, for orbital periods shorter than 4 days, 
the systems have circular orbits. Indeed, with such periods, the 
circularization of systems composed of O-type stars is generally 
already done or almost achieved. The mass ratio {Mp/Ms ) is uni- 
formly distiibuted between 1.0 a nd 10.0 and the mas s of the pri- 
mary is obtained from tables of Martins et al.l (12005) according 
to the spectral classification of the star. The longitudes of perias- 
tron are uniformly drawn between and 27r, and the orbital in- 
clinations are randomly drawn according to cos / e [-1; 1] from 
a uniform distribution. 

Table 4. Binary detection probability (%) for the time sampling 
associated to the different presumably single stars in our sample 
and for various period ranges (expressed in days). 



Stars Short Intermed. Long All 

[2-10] [10-365] [365 -3165] [2-3165] 

94.1 
86.6 
90.2 
87.8 
94.5 
78.8 
91.2 
84.0 
83.9 
86.0 
75.6 
92.4 
89.5 
86.4 
88.0 



The results of these simulations are listed in Table |4] This ta- 
ble gives the percentage of systems that would be detected with 
the same temporal sampling as our survey if their orbital pe- 
riod was short, intermediate, long or covering a timescale going 
from 2 to more than 3000 days, respectively. This reveals that 



at least 72% of the systems with periods smaller than one year 
should have been detected. However, between 1 1 and 64% of the 
systems with an orbital period larger than one year could have 
been missed. These results show that the detection efficiency of 
our survey is rather good, especially for systems with periods 
shorter than one year, but that we may also have missed some 
of the long-period binary systems. These results thus emphasize 
that the short-period systems are easier to detect than the long- 
period systems with this strategy of observation. 

7.2. Binary fraction in Cygnus OB associations 

The spectroscopic analysis made on the Cygnus region includes 
O-type stars coming from different environments. Therefore, 
considering them altogether boils down to study a random sam- 
ple of O-stars. It is also possible to focus on each association 
or even on young open clusters belonging to these associations. 
However, in the latter case, we would deal with small number 
statistics. The binary fraction will thus be discussed both over 
the entire region and for each OB association. 

We observed nineteen stars or multiple systems contain- 
ing at least one O-type object. Among this sample, one star, 
HD 229234, was classified as SBl and three were repotted 
as SB2 systems, HD 193443, HD 194649 and HD 228989. 
Therefore, no doubt remains on the binary status of these obj ects. 
As a first step, we focus on each individual OB association by 
first including only our sample stars. In CygOBl, the minimal 
binary fraction is of 33% (3 stars out of 9). In Cyg OB3, we find 
no star in a binary system out of 4, i.e., 0% whilst in CygOB8 
and CygOB9, we obtain 0% (0 out of 3) and 33% (1 out of 
3), respectively. By putting together these numbers to achieve 
a more global view of our sample, we find a minimum binary 
fraction of about 21% among the nineteen sample stars (4 out of 
19 ). However, i f we ta ke into account all the O stars mentioned 
bv 'Hum phrevsl (Il978h . the minimum O-type star binary frac- 
tions become equal to 33% (4 out of 12, BD +364063 also de- 
tected as binary with P^ib - 4.8 days,|W iniams et al. 2009), 33% 
(3 out of 9, HP 19120 1, iBurkholder et al . 1997, HD 1 90918 
lHilI&Underhiii|[T99l and HD 226868, [Oies & Bolton [1981 
were detected as binaries with orbital periods of 8.3, 112.4, 5.6 
days, respectively), 0% (0 out of 4) and 14% (1 out of 7) for 
the CygOBl, CygOB3, CygOB8 and CygOB9 associations, 
respectively. Overall, we thus obtain a minimal binary fraction 



HD 193514 


99.8 


96.0 


84.1 


HD 193595 


99.7 


86.6 


64.5 


HD 193682 


99.8 


92.4 


73.6 


HD 194094 


99.2 


87.0 


68.8 


HD 194280 


99.9 


96.0 


85.1 


HD 228841 


99.0 


81.4 


44.3 


HD 190864 


99.9 


94.5 


76.0 


HD 227018 


99.8 


82.9 


57.1 


HD 227245 


99.5 


68.3 


59.5 


HD 227757 


99.5 


82.4 


63.3 


HD 191423 


98.8 


72.3 


36.5 


HD 191978 


99.8 


94.1 


79.5 


HD 193117 


99.8 


91.0 


72.1 


HD 194334 


99.7 


91.3 


63.3 


HD 195213 


99.7 


89.2 


67.9 
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Fig. 26. Distributions of orbital parameters (period, eccentricity, mass ratio) of binary systems located in young open clusters or OB 
associations. 



of 25% (8 stars out of 32). Of course, this value was deter- 
mined by focusing only on the CygOBl, CygOB3, CygOBS 
and Cyg OB9 associations and is therefore not representative of 
the entire Cygnus region. Therefore, we include in this discus- 
si on the results of the lar g e spec troscopic/photometric survey 
of iKiminki & Kobulnickvl (20121 and the subsequent papers). 
These authors have indeed analysed 114 stars belonging to the 
CygOB2 association. This sample takes into account massive 
objects for which the primary stars are classified between B 3 
and 05. They found a hard minimum binary fraction of 21% 
(24 out of 1 14). Although our sample is clearly smaller, the min- 
imum binary fraction appears to be similar in most of the OB 
associations in the Cygnus complex. Finally, putting all these re- 
sults together, we reach a general minimum binary fraction of 
22% (32 out of 146) in the Cygnus region. 

7.3. Massive binaries in a wider context 

To bring additional constraints to this work in a larger frame, 
we have selected the SBl and SB2 systems (including at 
least one O star) detected in various spectroscopic multiplic- 
ity st u dies dSanae t al. 2008, 2009, 20111 iRa uw & De Becker 
2004"; 'De Becker et al. 2006: .Hillwig et alj " l2006; ,Rauw et aL 



2009; Mahy et al. 200^ as well as in numerous pa pers (see 
e.g., "spectroscopic" references in Mason et al. 1998) including 
studies devoted to the CygOB2 as sociation (Rauw et al. 1999: 
iDe Becker etall l2004t iNaze et al.1 120081: Kiminki eLalj lloTZ 
and the subsequent papers). Finally, we complete this dataset by 
adding the orbital information concerning the O-type st ars pro- 
vided by the 9* spectroscopic binary catalogue (IPourbaix et al.l 
[2004). Orbital solutions are then known for more than 140 sys- 
tems (SB 10 and SB2) which are mainly located in young open 
clusters or OB associations. We present, in Fig. |26] the corre- 
sponding P vs. e, P vs. q (= Mp/Ms) and q vs. e diagrams. In 
this figure, the filled symbols represent the SB2s, the open ones 
the SB 1 systems, and the squares give the parameters quoted in 
the 9* spectroscopic binary catalogue. In order not to affect the 
readibility of these diagrams and because some are missing, we 
do not include the error-bars on these measurements. 

From the period-eccentricity diagram (left panel of Fig. |26| ), 
we see that shorter period systems are dominantly character- 
ized by lower eccentricities, suggesting a lack of highly eccen- 
tric short period systems as well as almost circular long period 
systems. The majority of the datapoints seems to show a trend 
between the period and the eccentricity. However, the ques- 



^ SBl systems have not been included in the P vs. q and q vs. e plots 
because of the difficulty to infer reliable mass ratios in such systems. 



tion of the complete ness of the data has to be asked. Indeed, 
ISana & EvansI (1201 1*) mentioned that at least half the known and 
suspected spectroscopic binaries lack a reliable orbital solution. 
Among the systems that lack a reliable solution, the majority has 
long orbital periods. The apparent linear trend between eccen- 
tricity and orbital period could be affected by two competing ob- 
servational biases. Indeed, highly eccentric SB2 systems display 
large RV separations at periastron which are easier seen than the 
lower RV excursions of circular systems. However the duration 
over which these large separations are observed is rather short 
and could easily be missed. 

We see from the Pws.q plot (middle panel of Fig.l26]l that the 
majority of detected systems are those with a period smaller than 
30 days and with a mass ratio between 1 and 3. This plot does not 
agree with a uniform di stribution of the mass ratios in the range 
I < q < 5 as lSana & Ev ans (2011) suggested it. However, these 
authors also reported a decrease of the number of systems with 
q > 1.7 as it can be see in Fig.|26] These results clearly show the 
limitations of spectroscopy. Indeed, it is hardly possible to detect 
a secondary signature for systems with large mass differences 
without having high-quality high signal-to-noise data. 

Finally, the q vs. e plot (right panel of Fig. l26T l shows that 
a larger amount of systems with small eccentricities have mass 
ratios between 1 and 3 or said differently that only systems 
with almost similar components display a small eccentricity. 
However, the detectability of systems with a high mass ratio or 
a high eccentricity is difficult. Our sample could thus be biased. 
The decreasing range of high eccentricities for higher mass ra- 
tios, observed in the right panel of Fig. |26]could thus likely be 
stemmed from decreasing completeness, making this conclusion 
less reliable. Indeed, the tim e required for the circularization, 
given by Hu rlev et al.l (l2002h . is dependent on the mass ratio 
(q - (Mp/Ms)) squared. Therefore, when q is very high, the the- 
ory predicts that the time to circularize the system will be larger, 
which is not the case here. 

Although this sample remains small in comparison to the 
overall population of binaries among the O-type stars, these 
three diagrams provide a first approximation of the general dis- 
tribution of orbital parameters. It would be interesting however 
to include in these diagrams binaries which would be detected 
by adaptive optics, speckle interferometry or interferometry to 
cover larger orbital parameter distributions. 



8. Conclusion 

We revisited the binary status of nineteen O-type stars located in 
different OB associations of the Cygnus region. We confirm the 
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binarity of four objects, three SB2s: HD 193443, HD 194649, 
HD 228989 and one SBl: HD 229234. All these systems have 
short-term orbital periods, less than 10 days. We also derived 
for the first time the orbital parameters for the three SB2 sys- 
tems. The apparent lack of intermediate and long period sys- 
tems in these OB associations contrasts with the case of the 
NGC 2244 young open cluster where only one longer period sys- 
tem HD 46149 has been identified (HD 465 73 being located i n 
MonOB2 association and not in NGC 2244. iMahv et alj[2009[) . 
However, this difference seems linked to the sample. Indeed, 
among the four stars belonging to Cygnus OB associations 
which were already known as binaries and thus not included 
in our sampl e, HD 190918 is also a long period binary (PqA = 
112.4 days, iHill & Underiiililll995h . This result therefore does 
not allow us to point to different conditions for massive star for- 
mation. 

The stars contained in our sample were chosen because of 
their brightness and because, for most of them, their binary sta- 
tus has not yet been established. These nineteen stars thus con- 
stitute a random sample of O-type stars. Our strategy of ob- 
servations allowed us to reach a binary detection rate close to 
90% for periods up to 3165 days and to determine a minimum 
spectroscopic binary fraction in this sample of 21%. By includ- 
ing the stars of Humphreys' catalogue (1978) which belong 
to the OB associations studied in the present paper but which 
were not analysed because their binarity was already known, we 
reach a minimum spectroscopic binary fraction of 25%. When 
we consider all these associations separately, we obtain 33% of 
O-type stars in binary systems in CygOBl, 33% in CygOB3, 
0% in Cyg OB8 and 14% in Cyg 0B9. All these values can be 
completed by the minimum b inary fraction of 21% quoted by 
iKiminki & KobuMckrvl (20121) in CygOB2. Finally, when we 
take all these binary fractions together, we reach a binary frac- 
tion for the Cygnus region of about 22% (32 stars out of 146). 

In addition to the analysis of the RV variations, several ob- 
jects show significant variations of their line profiles. These vari- 
ations are mainly observed for the He 1 447 1 and He ii 4686 lines. 
Since not all the lines are affected, this implies that these vari- 
ations are probably not due to binarity but rather result from 
phenomena intrinsic to stellar winds or perhaps even from non- 
radial pulsations. The spectral classifications of the stars derived 
in the present paper and the results of the TVS analysis are sum- 
marized in Table |5] 

The results presented here only focus on a small sample of 
O-type stars in the Cygnus complex. This region represents a 
large panel of O-type stars which can provide numerous clues 
on the formation scenarios of these objects. In paper II, we will 
continue the analysis of these stars by determining with a model 
atmosphere code their individual parameters, their age, and their 
chemical enrichment to better constrain the properties of these 
objects. 
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Table 5. Summary of the multiplicity of stars and of their line profile variations 



Stars 


Sp. type 


He I 4471 


TVS variations 
He II 4542 Hen 4686 


He 1 4713 


HD 193443 


09 + 09.5 


Y 


Y 


Y 


Y 


HD 193514 


7-0 7.5 Ill(f) 


Y 


Y 


Y 


N 


HD 193595 


07V 


Y 


N 


N 


N 


HD 193682 


5111(f) 


N 


N 


Y 


N 


HD 194094 


8in 


- 






- 


HD 194280 


09.7 


Y 


N 


N 


N 


HD 228841 


07n 


N 


N 


N 


N 


HD 228989 


8.5V + 09.7V 


Y 


Y 


Y 


Y 


HD 229234 


9III + ... 


Y 


Y 


Y 


Y 


HD 190864 


6.5111(f) 


Y 


N 


Y 


N 


HD 227018 


6.5V((f)) 


Y 


N 


N 


N 


HD 227245 


7V-lII((f)) 










HD 227757 


09V 










HD 191423 


ON9inn 


N 


N 


N 


N 


HD 191978 


08in 


Y 


N 


Y 


N 


HD 193117 


09in 


Y 


N 


N 


N 


HD 194334 


07-07.5 Ill(f) 


Y 


N 


Y 


N 


HD 194649 


6111(f) + 08V 


Y 


Y 


Y 


Y 


HD 195213 


07111(f) 


Y 


N 


Y 


N 



Notes. Notes: "Y" means that significant TVS variations in the spectral lines are observed whilst "N" reports the absence of variations. 
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